Orientation Dependence of the Intrinsic Anomalous Hall Effect in hep Cobalt 
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We carry out first-principles calculations of the dependence of the intrinsic anomalous Hall con- 
ductivity vector it" of hep Co on the magnetization direction m. The magnitude of cr" decreases 
smoothly from 481 S/cm to 116 S/cm as the magnetization is tilted from the c-axis to the a6-plane. 
This factor-of-four reduction compares well with measurements on single crystals, while the angular 
average (cr" ■ rh) = 226 S/cm is in excellent agreement with the value of 205 S/cm measured in 
poly crystalline films. The strong anisotropy of ct" is a consequence of spin-orbit induced changes in 
the electron states near the Fermi level as the magnetization is rotated. 



PACS numbers: 75.30.Gw,75.30.-m,72.15.Gd,78.20.Ls 

The Hall effect in ferromagnets has two contributions. 
The ordinary Lorentz- force part, which is proportional 
to the magnetic field for weak fields, and an "anoma- 
lous" part, which depends on the magnetization. A the- 
ory of the anomalous Hall effect (AHE) was put forth 
by Karplus and Luttinger [1], who showed that a Hall 
current perpendicular to the electric field and odd under 
magnetization reversal is established in a ferromagnetic 
crystal as a result of the spin-orbit interaction (SOI). 

The link between the AHE and the SOI became con- 
sensual, but a long debate ensued on whether the rele- 
vant SOI is associated with the crystalline potential (in- 
trinsic) or with impurity atoms (extrinsic). The asym- 
metric impurity scattering of the spin-polarized charge 
carriers (skew-scattering) leads to a linear dependence of 
the anomalous Hall resistivity — Ey/jx on the longi- 
tudinal resistivity p^x] an additional scattering process, 
side-jump, yields the same scaling Py^ cx p^.^ as the in- 
trinsic Karplus-Luttinger mechanism. For a review see 
Ref. 0. 

Unlike the extrinsic contributions, which depend on 
the details of the impurity potential, the Karplus- 
Luttinger anomalous Hall conductivity (AHC) can be 
evaluated from the band structure of the crystal 0, |3| : 
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where /„k is the Fermi-Dirac distribution function and 
n„^y (k) is the Berry curvature tensor of each cell- 
periodic spinor Bloch state kt„k)- First-principles cal- 
culations for Fe, Co, Ni [E i, H, SrRuOg @ and 
MnsGesQ have consistently found good agreement with 
room-temperature experiments, establishing the impor- 
tance of the intrinsic mechanism. 

Recent experiments have focused on isolating the dif- 
ferent contributions to the AHE 0,11, 0]. Skew-scattering 
can be separated from the other two terms by fitting the 
measured anomalous Hall resistivity to the form 



where b — a^y 
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It''^ . The coefficients a (skew-scattering) 
and b (intrinsic plus side-jump) can be read off a plot of 
Pyxl Pxx versus p^x, where p^x is varied through doping 
or temperature changes. 

In this Letter, we present a detailed first-principles 
study of the orientation dependence of the intrinsic AHE. 
Anisotropy in the AHE has been measured in single crys- 
tals of the ferromagnetic elements (hep Co [10| , hep Gd 
11 1, fee Ni [l2|, and bcc Fe flslV and more recently in 
ferromagnetic compounds [l4L llsL [l6l | . On the theoret- 
ical side there has been little progress beyond the ba- 
sic phenomenological description. To our knowledge, the 
only attempt at a microscopic model has been the tight- 
binding study of Ref. |l4|. Because the AHC is very sensi- 
tive to fine details in the band structure [1,|3,[11j a quanti- 
tative ah initio theory of anisotropy is highly desirable. It 
is also not obvious that the phenomen olog ical description 
of magnetocrystalline anisotropy IT], ll8| applies to the 
AHC given by Eq. ([T|). The Berry curvature undergoes 
strong and rapid variations in A:-space, with sharp peaks 
and valleys from avoided crossings near the Fermi level 
0, 0, 0| • It has been argued that such behavior cannot be 
described perturbatively Q, and that it often gives rise 
to a complex or even irregular behavior of the AHC as 
a function of exchange splitting and Fermi level position 
. This raises the possibility that the orientation depen- 
dence of (T°^- may also not be smooth. We find instead 
that in hep Co it is strong but remarkably smooth, and 
can be described by a phenomenological power-series ex- 
pansion. The calculated anisotropy accounts for the ex- 
perimental observations in both single crystals (angular 
dependence) and polycrystalline films (angular average). 

We begin by reviewing the phenomenology 17, [l^ . 
Electrical conduction in ferromagnets is described by 
a magnetization-dependent conductivity tensor: Ji = 
aijCM^Ej, where i,j are Cartesian indices. The Onsager 
relation CTij(M) = CTji(— M) implies that the symmetric 
{afj ) and antisymmtric {afj ) parts of aij are respectively 
even and odd functions of the magnetization M. The cur- 
rent density J is therefore comprised of an even (Ohmic) 
part J'* and an odd (Hall) part J°. The Hall current 
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FIG. 1: (Color online.) Evolution of the components of the 
anomalous Hall conductivity parallel (um) and perpendicular 
{ag) to the magnetization [Eq. ([Sjl] as M is tilted by 6 from 
the c-axis towards the a-axis. The solid lines are fits to the 
first-principles data, as described in the text. The left and 
right insets shows respectively (Tm(7r/2, ip), 



and (7^(6', 0). 



reads J" = E x cr", 



where cr£ 



{l/2)e,jkaij. 



IS per- 



pendicular to E but not necessarily to M, since cr° and 
M may not be coUinear. This is the situation in sin- 
gle crystals, where cr" || M only when M points along 
certain high-symmetry directions. 

To characterize the anisotropy we write cr" as 



(3) 



where 6 and (p are respectively the polar and azimuthal 
angles of M = Mm. In macroscopically isotropic sys- 
tems such as polycrystals cr" || M, i.e., (jq = a^p = 0. The 
orientation dependence in hep single crystals is given to 
third order in a spherical-harmonic expansion by [l9| 



O-f = CiiFii + 031^31 
erf = CiiFi,-i + 031^3, 
erf = Cioi^lO + C3oy30 



(4) 



where Yim{d,(p) are real spherical harmonics. Because 
erf is independent of Lp, while cr" and ctj have respec- 
tively cosine and sine dependences, tr" and M share the 
same azimuthal angle, and their polar-angle mismatch is 



independent of Lp. Thus cr„ 



fTm,9iO) and = 0. 



We have carried out fuUy-relativistic band-strucutre 
calculations for hep and fee Co at the experimental lat- 
tice constants of 4.74 and 6.68 bohr respectively, using 
the PWSCF code [2^ . The pseudopotential was generated 
using similar parameters as in Ref . |3- The plane- wave ba- 
sis cutoff was set at 140 Ry, and the exchange-correlation 
functional, fc-point sampling, and Fermi smearing used in 
the self-consistent calculation are the same as in Ref. Ijy. 
The integral in Eq. ((1) was evaluated using a Wannier- 
interpolation scheme ^ to sample efficiently the Bril- 
louin zone over a 125 x 125 x 125 uniform /c-point mesh 



TABLE I: Anomalous Hall conductivity \a°'\ = am in S/cm 
for selected high-symmetry orientations of the magnetization 
in hep and fee Co. The AHC of poly crystalline samples is 
calculated as an orientational average (see text). 



Co 


Orientation 


Calc. 


Expt. 


hep 


c-axis 


481 


683" 




a?)- plane 


116 


232" 




Polycrystal 


226 


205*, 275" 


fee 


[001] 


249 






[110] 


218 






[111] 


234 





"Reference [10, 
''Reference ^ 
'^Reference 23 



(200 X 200 X 200 for fee Co), with a 5 x 5 x 5 adaptively 
refined mesh around the points where the magnitude of 
the Berry curvature exceeded 10 A^. 

The AHC of hep Co was calculated for several ori- 
entations of the cell-averaged magnetization [2l| in the 
ac-plane {ip = 0). The tilting angle 9 was increased 
from (M |) c-axis) to tt/2 (M ]] a-axis) in steps of 
7r/32, and for each step the vector cr°-{9^p) was calcu- 
lated. Fig. [1] contains the numerical results: (Tm(^,0) 
and 'Jq{9, 0) are shown in the main panels, while the in- 
sets contain additional data which confirms the absence 
of (or very weak) basal-plane anisotropy. The vectors cr" 
and M start out parallel, but as M begins to tilt away 
from the c-axis <t° lags behind {ag < 0), and they be- 
come parallel again upon reaching the basal plane. The 
AHC is strongly anisotropic, decreasing in magnitude by 
a factor of 481/116 ~ 4.1 between 6* = and 6* = 7r/2. 
This is in reasonable agreement with the ratio of 2.93 
measured in single crystals at 290 K [13]. While strong, 
the angular dependence of cr° is smooth, and can be de- 
scribed by Eq. A least-squares fitting to the data 
yields, in S/cm, cio — 951.5, cn = —204.1, C30 = 1.2, 
and C31 — 38.4, producing the sohd-hne curves in Fig. [T] 

We now turn to the cornparison with the measurements 
on polycrystalline films [8|. The films were magnetized 
along the growth direction by an applied field; assum- 
ing randomly oriented crystallites, each with a bulk-like 
Hall current density J"' {9, p), the net Hall current density 
in the films can be estimated by performing an orienta- 
tional average: (J°) = E x (cr°). Because cr° displays 
azimuthal isotropy, it suffices to average Eq. ([3]) over 9 
for fixed p. The average of (Jg vanishes (see Fig. [T]), re- 
sulting in an isotropic AHC (cr") = {am)rh of magnitude 

^poiy = lo^^ a^{9)smede = 226 S/cm. This should be 
compared with the value h — 205 S / cm obtained in Ref. d 
by fitting the experimental data to Eq. ((2|). The exper- 
iment does not discriminate between the intrinsic and 
side-jump components of 6, but the close agreement with 
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FIG. 2: (Color online.) Upper panel: MCD spectrum for two 
magnetization directions. Lower panel: Cumulative contribu- 
tion to the AHC from the spectrum above energy huj. The 
circles denote the AHC calculated directly from Eq. ((TJ. In- 
set: Cumulative contribution to the self-rotation part of the 
orbital magnetization from the spectrum below hu:. 



FIG. 3: (Color online.) Upper panel: Energy bands close 
to the Fermi level for M || c-ajds (0 = 0) and M || a-axis 
(9 = 7r/2). Middle and lower panels: fc-space Berry curvature 
nk(6') for e = 0, e = tt/2, and e = 7r/4. The full height of the 
resonance peaks near the L-point is of the order of 10'^ A^. 



the calculated o-^o\y reinforces the conclusion 8] that the 
former dominates. Table U summarizes the comparison 
between our calculations and experiments. 

Next we discuss the origin of the strong anisotropy. 
The AHE of uniaxial crystals is anisotropic to first order 
in an expansion in powers of the magnetization To^ , while 
in cubic crystals anisotropy appears only in third-order 
[l^ . and is expected to be much weaker. For example, 
the AHC of fee Co changes by less than 10% as a function 
of the magnetization direction (Table Perhaps more 
surprising is the fact that the AHE in hep Co appears to 
be considerably more anisotropic than both the magneto- 
optical spectrum and the orbital magnetization 24, 
25| . This is intriguing because the three phenomena are 
related by linear sum rules (26j . and hence anisotropy 
appears at the same order. 

The sum rules read {uj-Hma-")^ = (7r/2)cr'^(a; = 0) 
and {lm(T°-)^ = (7rec/?i)Mg^, where {f)uj = f{uj)duj. 
^"•{uj = 0) is the dc AHC; at finite frequencies cr" ac- 
quires an imaginary part which describes the differential 
absorption of right and left circularly-polarized light, or 
magnetic circular dichroism (MCD). The first sum rule 
expresses the AHC in terms of the first inverse moment 
of the MCD spectrum. The second relates Mg^, the 
"gauge-invariant self-rotation" part of the orbital mag- 
netization [26|, to the zero-th spectral moment. 

The absorptive part of (Tm(w) is plotted in the upper 
panel of Fig. [5] for = 0, 7r/2. The lower panel shows 
A^'^{uj) = I ^Imcr,„(w')rfw'. For t^max oo (in 

practice we use o^max = 7 eV) A^{0) — (Tm(w — 0), 
so that A^{uj > 0) is the cumulative contribution to 
the AHC from optical transitions above uj. While for ei- 
ther orientation there are sizeable contributions to the 
AHC up to w ~ 3.5 eV, its orientation dependence is 



concentrated below 0.3 eV. At these low frequencies the 
MCD spectrum changes sign between 9 = and 9 = tt/2. 
This difference gets magnified in the AHC via the inverse- 
frequency weight factor, producing the bifurcation below 
0.3 eV of the two A^^^co) curves. All frequencies are 
equally weighted in the orbital moment sum rule, which 
as a result is more isotropic. This is seen in the inset, 
where we plot Bf^{u;) = lma^{u;')duj' , the cu- 

mulative contribution below lo to the gauge-invariant self- 
rotation per atom. 

The origin of the low-frequency anisotropy can be seen 
in Fig. O The upper panel displays the energy bands 
for the two magnetization directions. Rotating M from 
the c-axis to the a-axis in the presence of the SOI turns 
various band crossings into avoided crossings and vice- 
versa. When this occurs close to the Fermi level the 
Berry curvature along M can flip sign in the process, 
while retaining a large magnitude. This is what happens 
near the L-point, as seen in the middle and lower panels, 
where we plot — fik • A for the two orientations (fik.fc = 
(l/2)eyfc J2n fnk^nk,ij IS the total Berry curvature at k). 
The sensitivity of the AHC to changes in the electron 
states near Ep may also be understood from the fact 
that cr° can be recast as a Fermi surface integral [El, , 
whereas the orbital magnetization truly depends on all 
occupied states. 

How can the spiky behavior of Jlk be reconciled with 
the smooth angular dependence displayed by a"' (9) in 
Fig. [1]? According to the phenomenological expansion 
dH), erf (X Mi to leading order (i = x,y,z). This will 
be the case for the AHC given by Eq. ([1]) provided that 
ilk.i oc Mi at each k. This proportionality holds reason- 
ably well even around strong resonance peaks, judging 
from the comparison in Fig. [3] between fik calculated at 
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anisotropy of the skew-scattering contribution [16'|. 
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FIG. 4: Evolution of cr" as the magnetization is rotated in the 
ac plane, plotted as ct" (cos 61) and (t" (sin ^). Inset: Anisotropy 
ratio a'^{0 — 0)/a^{6 = 7r/2) versus the reduced magnetiza- 
tion, according to the spin-fluctuation model. The dotted line 
denotes the approximate location of the hcp^fcc transition. 



9 = 7r/4 and at 6* = 0,7r/2. 

We end with a discussion of temperature effects. The 
AHC hardly changes as the Fermi-smearing tempera- 
ture in Eq. IJ) is varied from K to 300 K (Fig. [I]). 
This agrees with the constancy of the coefficient b of 
polycrystalline films in the range 78-350 K [8]. The 
angular dependence of the AHC can also give rise to 
a temperature dependence, via long-wavelength ther- 
mal fluctuations in rh Q. According to this model, if 
erf (T — 0) changes linearly with Mi upon rotating M 
then (T'^(T) = [M(r)/M(0)]cr'^(r = 0), which seems 
to explain the temperature dependence in MnsGea Q- 
Fig. |4] shows that in hep Co cr" depends linearly on , 
while the a^^M^) curve is significantly nonlinear (note 
that |c3i/cii| » |c3o/cio|). As a result, the a-axis AHC 
should decrease with T faster than M{T), producing an 
increase with temperature of the ratio a'^{9 — Q)/a^{9 — 
tt/2). An estimate of the magnitude of this effect can be 
obtained using the "Z(Z + l)/2£0wer law" for the 

coefficients cim{T) (see Ref. 'l9 for details). The result, 
shown in the inset of Fig.dl is a 17% increase between K 
and Tc- This is a sufficiently large effect that is should be 
observable in principle, but in practice it is preempted by 
the phase transformation into the fee structure at 695 K, 
wen below Tc = 1400 K 0. 

In summary, we have shown by means of first-principles 
calculations that the intrinsic mechanism for the AHE 
describes quantitatively the observed strong angular de- 
pendence in hep Co single crystals. The key role of 
near-degeneracies across the Fermi level was elucidated, 
and the AHE of polycrystalline Co films was reproduced 
by averaging the single-crystal Hall conductivity over all 
magnetization directions. Further experimental and the- 
oretical studies of the orientation dependence of the AHE 
are needed. For example, very little is known about the 
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